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Abstract —Massive multiple-input multiple-output (MIMO) is a 
promising technology aiming at achieving high spectral efficiency 
by deploying a large number of base station (BS) antennas using 
coherent combining. Channel aging due to user mobility is a 
significant degrading factor of such systems. In addition, cost 
efficiency of massive MIMO is a prerecuisite for their deployment, 
that leads to low cost antenna elements inducing high phase 
noise. Since phase is time-dependent, it contributes to channel 
aging. For this reason, we present a novel joint channel-phase 
noise model, that enables us to study the downlink of massive 
MIMO with maximum ratio transmission (MRT) precoder under 
these conditions by means of the deterministic equivalent of the 
achievable sum-rate. Among the noteworthy outcomes is that 
the degradation due to user mobility dominates over the effect of 
phase noise. Nevertheless, we demonstrate that the joint effects of 
phase noise and user mobility do not degrade the power scaling 
law I/s/m {M is the number of BS antennas), as has been 
established in massive MIMO systems with imperfect channel 
state information. 

Index Terms —Channel estimation, Doppler shift, phase noise, 
massive MIMO, linear precoding. 

I. Introduction 

The large multiple-input multiple-output (MIMO) paradigm, 
known also as massive MIMO lU, is a promising technology 
for achieving higher throughput per area, which is necessary 
for next generation systems 121 ■ In particular, this concept 
employs multi-user MIMO technniques, where the base station 
(BS) includes a large antenna array that results to the vanish¬ 
ment of both inter-user interference and thermal noise as the 
number of antennas increases due to channel orthogonality. 

Originally, most studies assume perfect channel state infor¬ 
mation (CSI); however, this is a strong assumption, being 
untenable in practice. Given that the performance depends 
on the quality of CSI at the transmitter, it is of paramount 
importance to investigate realistic scenarios with imperfect 
CSI, acquired during the uplink training phase in time-division- 
duplex (TDD) architectures fTl. 121. 

One of the main sources contributing to imperfect knowl¬ 
edge of CSI is channel aging, present in time-varying channels, 
that describes the mismatch between the estimated and current 
channels due to the relative movement between users and 
BS antennas. Despite its crucial role, it has not been studied 
in depth. Specifically, the uplink and downlink signal-to- 
interference-and-noise-ratios (SINRs) have been derived by 
using a deterministic equivalent analysis after applying a 
maximum ratio combining (MRC) decoder and a maximum 
ratio transmission (MRT) precoder, respectively a. Furter in¬ 
vestigation and comparison took place in ||5], ||6|, where more 
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sophisticated linear techniques were applied. Moreover, the 
optimal linear receiver for cellular systems has been derived in 
171 by exploiting the correlation between the channel estimates 
and the interference from other cells, and in JS], the uplink 
analysis of a cellular network with zero-forcing (ZF) receivers 
for any finite as well as infinite number of BS antennas has 
been provided. Neverteless, tight lower bounds and the power 
scaling law have been studied for massive MIMO with channel 
aging in Q. 

In practice, the application of massive MIMO is meaningful, 
if both hardware cost and power consumption are kept low. 
In such case, hardware impairments become a significant de¬ 
grading factor leading to imperfect CSI. Among the hardware 
imperfections, phase noise ITOl . randomly phase shifting the 
output of the local oscillator (LO), accumulates over time. Un¬ 
fortunately, it has been scarcely studied in the case of massive 
MIMO systems, and mostly for uplink scenarios. Specifically, 
the phase noise has been taken into account for single carrier 
uplink massive MIMO by considering time-reversal MRC and 
ZF in ifm and ifTTl . respectively. Furthermore, in llT3ll . the 
achievable uplink user rate with MRC has been obtained by 
considering a general model with additive and multiplicative 
(phase noise) distortions, and it was shown that the phase 
noise can be alleviated in the large antenna limit. Especially, 
regarding the downlink, an effort to address the effect of phase 
noise, being of great importance in large MIMO systems, has 
taken place only in fMl . 

Given that phase noise, being time dependent, contributes 
to channel aging, this work provides a novel model incor¬ 
porating both effects of user mobility and phase noise. In 
fact, we present an insightful joint channel-phase noise model. 
In addition, we provide the deterministic equivalent of the 
downlink SINR with MRT. Furthermore, the power scaling 
law has been studied in the case of user mobility combined 
with phase (general channel aging conditions). Numerical 
results reveal the impact of both effects in massive MIMO 
for various mobility conditions and topologies of LOs. Among 
the interesting observations, we mention that the effect of user 
mobility is more severe than phase noise, and the massive 
MIMO technique keeps being favorable even in general time- 
varying conditions. 

II. System Model 

We consider a single-cell system with a BS serving K 
single-antenna non-cooperative users. In particular, the BS is 
deployed with an array of M antennas. Our study focuses 
on large scale networks, where both the number of antennas 
M and users K grow infinitely large, while keeping a finite 
ratio P, i.e., M,K -A oo with K/M = P fixed such that 
limsupjv^AT/M < oo. 
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A realistic quasi-static fading model is taken into account, 
where the M x 1 channel response vector between the BS 
and user k during the nth symbol, is denoted by hk^n S C^. 
Especially, „ G exhibits flat fading, and it is assumed 
constant during the current symbol, but it may vary slowly 
from symbol to symbol. Moreover, common effects such as 
path loss dependent on the distance of the users from the BS, 
and the spatial correlation due to limited antenna spacing, are 
considered by expressing the channel between as 

hk,n = ( 1 ) 


where Rt = E 




G C- 


MxM 


a deterministic 


Hermitian-symmetric positive definite matrix representing the 
aforementioned effects. The independence of R^ by n is a 
reasonable assumption, since effects such as spatial correlation 
and lognormal shadowing vary with time in a much slower 
pace than the coherence time a, a, Q. In addition, „ G 
is an uncorrelated fast fading Gaussian channel vector 
drawn as a realization from zero-mean circularly symmetric 
complex Gaussian distribution, i.e. ~ CA/'(0 ,Im). 


A. Phase Noise Model 

In practice, phase noise, induced during the up-conversion 
of the baseband signal to bandpass and vise-versa, impairs 
both the transmitter and the receiver, and it appears as a 
multiplicative factor to the channel vector. The cause is a 
distortion in phase due to the random phase drift of the 
generated carrier by the corresponding LO. 

Our analysis regarding the nth time slot, takes into account 
for both synchronous and non-synchronous operation at the 
receiver. The latter setting assumes independent phase noise 
processes (j)m,n,'m = with being the phase 

noise process at the mth antenna. Note that the phase noise 
processes are considered as mutually independent, since each 
antenna has its own oscillator, i.e., separate independent LOs 
(SLOs) at each antenna. This setting can degenerate to identi¬ 
cal phase noise processes, if the LOs are assumed identical 
(ILOs). In the synchronous scenario, all the BS antennas 
are connected with a common LO (CLO), i.e., there is only 
one phase noise process Lurthermore, we denote (pk,n, 
k = 1,... ,K the phase noise process at the single-antenna 
user k. 

Phase noise during the nth symbol can be described by 
a discrete-time independent Wiener process , i.e., the phase 
noises at the LOs of the mth antenna of the BS and fcth user 
are modeled as 


0m,n — 0m,n —1 “t” Und ^k,n — Pk,n — 1 ; (^) 

where <5^™ ~ A/'(0,cr^^) and ~ A/’(0,cr^^). Note that 
af — in^fcCiTs, i = <j>m,Pk describes the phase noise 
increment variance with Tg, Ci, and /c being the symbol 
interval, a constant dependent on the oscillator, and the carrier 
frequency, respectively. 


The knowledge of the downlink channel is considered known 
due to the property of reciprocity. 

By assuming that the channel estimation takes place at time 
0 , we now derive the joint channel-phase noise linear minimum 
mean-square error detector (LMMSE) estimate of the effective 
channel o = ®fc,ohfe,o under phase noise, small-scale 
fading, and channel impairments such as path loss and spatial 

correlatiorQ. Note that 0fc,o = diag ..., | is 

the phase noise because of the BS and user k LOs at time 0 

(i) 

with 01'q = + Pk,o, i = l,...,M. 

Proposition 1: The LMMSE estimator of gfc,o, obtained 
during the training phase, is 

gM=(iM + ^R,') (3) 

where q is a noisy observation of the effective channel 
from user k to the BS, and pp = rp^ with p^ and r being 
the power per user in the uplink data transmission phase and 
the duration of the training sequence during the training phase, 
respectively. 

Proof: During the training phase T, all users transmit 
simultaneous mutually orthogonal pilot sequences consisting 
of length T symbols. Also, we assume that the effective chan¬ 
nel remains constant during this phase. In particular, the pilot 
sequences can be represented by fl' = [- 01 ; • • • ; 't^K] G 
with normalized, i.e., = I^. Given that the channel 

estimation takes place at time 0, the received signal at the BS 
is written as 


Yp,o =-v/^0fc,oHofI' + Zp_o, (4) 

where pp is the common average transmit power for all users, 
and Zp 0 G is spatially white additive Gaussian noise 

matrix at BS during this phase. By correlating the received 
signal with the training sequence of user k, the BS 

obtains 

y k,o ~ ®fc,oh/c,o “b ^=2p Q, (5) 

vPp 

where Zp_o = 'Zp^ipk ^ CAbjO, ct^Im). Since we seek to 
And the joint channel-phase noise estimate, we denote the 
effective channel at time 0 to be g^o = 0fe,ohfe,o- After 
applying the MMSE estimation method M, the effective 
channel estimated at the BS can be written as in Q. 

Employing the orthogonality principle, the current channel 
decomposes as 

gfe,o = gfe,o + gfc.o, ( 6 ) 

where g^o is distributed as CAb(0,Dfc) with = 

(Im + Rfc’ and gk,o ^ CJ\f{0, Rfc - D^) is the 

channel estimation error vector. Note that gk and g^ are 
statistically independent because they are uncorrelated and 
jointly Gaussian. ■ 


B. Channel Estimation 

In real systems obeying to the TDD architecture, the BS es¬ 
timates the uplink channel 0, 0, IISl. The estimation occurs 
by means of an uplink training phase during each transmission 
block, where the transmission of pilot symbols takes place. 


C. Channel Aging 

The relative movement between the users and the BS 
results to a Doppler shift, which makes the channel vary 

'Since the duration of the training phase is small, the phase noise can be 
considered constant because its innovation is negligible. 
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between when it is learnt by estimation, and when the estimate 
is applied for detection/precoding. The higher the velocity, 
the greater the performance loss because of worse channel 
estimation. The description of channel aging can be made by 
an autoregressive model of order 1, which relates the joint 
channel-phase noise process „ between the BS and the fcth 
user at time n (current time) to the estimated channel during 
the training phase. Specifically, we have 

Sk,n — ^nSk,0 “b (7) 

where, initially, A„ = Jo(27r^rsn)lM is supposed to model 
2-D isotropic scattering ifTTIPl gfc o is the channel obtained 
during the training phase, and is the uncorrelated 

channel error vector due to the channel variation modelled as 
a stationary Gaussian random process with i.i.d. entries and 
distribution CJ\f(0, — A„RfcA„). 

Since (|6]l allows to expresss the channel at time 0 by means 
of its estimate, the effective channel at time n can be written 
as 

S/c,n —-^nSk,0 “b ^k,n 

~-^nSk,0 “b ^k,n: (8) 

where ek,n = Angk,n + ek,n ^ CJ\f(0, Rk - A„DfeA„) and 
gfc 0 are mutually independent. As far as A„ is concerned, 
it is assumed known by the BS. Interestingly, it incorporates 
the effects of 2-D isotropic scattering, user mobility, and phase 
noise, as shown by Theorem[T] Consequently, useful outcomes 
can be extracted during the analysis in Section |III] 

Theorem 1: The estimated effective joint channel-phase 
noise channel vector at time n is expressed as 

gfeT* = A"gfe,o (9) 


Remark 1: Notably, phase noise induces an extra loss in 
the coherence due to its accumulation over time between the 
actual channel and the estimated channel coefficients. In other 
words, phase noise is another contribution to the channel aging 
phenomenon that imposes a further challenge to investigate the 
realistic potentials of large MIMO. 

Corollary 1: In case that all BS antennas are connected 
with the same oscillator (CLO) or all the oscillators 
are identical (ILO), A„ degenerates to a scalar = 

Jo(27r/£)Tsn)e 2 " or a scaled identity matrix arJ-M, 

respectively. 

D. Downlink Transmission 

In dowlink, the BS broadcasts data to its serving users with 
the same power pd- Exploiting that in TDD the downlink 
channel is given by the Hermitian transpose of the uplink 
channel, the received signal £ C by user k during the 
data transmission phase {n = t + 1,. .., T^) is 

yk,n — s/Pd^k,n^k,n^n “b (12) 

where s„ = a/AF„x„ denotes the transmit signal vec¬ 
tor by the BS with A, F„ £ and x„ = 

[a;i,n, X2,n,--- , XK,ny £ -- CA/'(0,Ik) being a 

normalization parameter, the linear precoding matrix, and the 
data symbol vector to its K serving users, respectively. Also, 
Zk,n ^ complex Gaussian noise at user k. 

The normalization parameter is obtained by constraining the 
transmit power per user to pd, i.e., E s„] = pd, as 

A = —n— -T- (13) 

E[itrF„FH] 

Since gk,n = ®k,niik,n, the received signal is written as 


with 

A„ = Jo(27r/z5r,n)e-—(10) 

r 

where = diag < e 2 e 2 

Proof: The total mean square error (MSE), which is 
actually given by means of the error covariance matrix 'Ek,n, 
is MSE = trEfc Specifically, 'Ek,n can be written as 



yk,n — ''J ^Pdgk,n® k,n^k,nXk,n 

+ \/^PdSk,n®k,n^i,nXi,n + Zk,n 

i^k 

= \/k,n^k,7i\ ^k,n~^^k,n 

desired signal noise 

+ ^ ^P<igk.n®k.n^i,'nXi,n 
i^k 


(14) 


Ffc,n — E [(gfc,n Ayigfc q) {gk,n -^ngkp) ] 

= IE|gfc,„P + A„E|gfe,opA,^j - 2E [Re {A„gfe,ogfc,„}] 

= Rfc + A„RfeA" - 2E [Re {A„gfc,ogfe,n}] 

2 

= Rfc -b A„RfeA" - 2Jo(27r/£,Tsn)e 5-"A„RfcA$„, 

( 11 ) 

where we have used that gk,n = ®fc,nhfc„ with &k,n 
and hfc „ being uncorrelated. In addition, we have denoted 

= diagje ^",...,e and E[hfc,oh^_„] = 

Jo(27r/DTsn)Rfc. An is determined by minimizing the MSE 
in Thus, if we differentiate with respect to A„, 

and equate the resulting expression to zero, we obtain the 
Hermitian matrix A„ as in (Uni). ■ 

^Normally, the BS antennas are deployed in a fixed space instead of 
considering distantly distributed small BSs constituting one large MIMO BS. 
In such case A„ is a scalar, since it is reasonably assumed that all antennas 
appear the same relative movement comparing to the user. 


Ry/Apd {Sk,n®k,n^k,nXk,n ^ [Sk,n®k,n^k,n\ Xk,n} ■ (15) 
We reasonably use the above technique because user k does 
not have access on instantaneous CSI, but he can be aware 
of ¥,[g'l Moreover, since the available CSI at time n 

is A„gfcp, as can be seen from (|8ll, the MRT precoder is 
Fn Aj^Gq. 

By taking into account that the input symbols are Gaussian, 
and by considering the worst case uncorrelated additive noise 
with zero mean as in na, we consider a SISO model, scaled 
by 1/M^, with desired signal power 


Sk,n — 


M2 


E 


Sfe,n®fe,nAngfc ,0 


and interference plus noise power at user k given by 


r ^ 

Ik,n = ^var 

A 


Sk,n®k,n^k,r, 


-b 


PdM2 


i^k 


Sk,n^k,n^i^r 


(16) 


(17) 

















4 


where we have denoted ©fc^„ = A©^ Note that we 
have assumed that Qgk,n, since gk,n is circularly 

symmetric. 

The mutual information between the detected symbol and 
the transmitted symbol is lower bounded by the achievable 
rate per user given by 

^ Tc — T 

= (1 + 7fc,«), (18) 

n—1 

o 

where ^k,n = is the instantaneous downlink SINR at time 
n, and the expectation operates over all channel realizations. 


III. Performance Analysis 
The investigation of the downlink SINR with MRT precoder 
accounting for the effects of imperfect and delayed CSIT due 
to user mobility and phase noise is the subject of this section. 
Specifically, we derive the deterministic equivalents of the 
SINR, being tight approximations even for moderate system 
dimensions. Simulations in Section HV] corroborate our results 
and their tightness. 

The deterministic equivalent of the SINR 7 ^ „ is such that 
Ik n — Ik n > (fl while the deterministic rate of user k 

’ ’ M—¥00 

is obtained by the dominated convergence and the continuous 
mapping theorem ED as 


Rk 


1 ^ ^ 

Y log2(l +7fc.») 

^ n—1 


a.s. 


^ 0 . 


N—^oo 


(19) 


The deterministic equivalent downlink achievable rate of fcth 
user with MRT follows. 

Theorem 2: The deterministic equivalent of the downlink 
SINR with MRT precoding, accounting for imperfect CSI due 
to general impairments, and delayed CSI due to phase noise 
and user mobility is given by 



4 = ^trA^Dfc, (5^ = j^tr A^Dfc (Rfc - A^Dfc), and 
= ]gtr A^D.Rfe. 

Proof: The proof of Theorem|2]is given in AppendixiAl ■ 
Proposition 2: For the identical SLOs or CLO setting, the 
downlink achievable SINR of user k with MRT precoding at 
time n subject to imperfect CSI, delayed CSI due to phase 
noise and user mobility, and collocated BS antennas becomes 

= ^J^^Jo(27r/Drsn)e-"("'»=+"^)"[R2]™„, (21) 

where the transmit uplink and downlink powers are scaled 
proportionally to i.e., andpd = E^/M’^ 

for fixed Eu^ and E^, and q > 0. 


^Note that 


denotes almost sure convergence, and a-n 


expresses the equivalence relation an — bn 
t\vo infinite sequences. 


-4 0 with an and bn being 


Proof: Let us first substitute the MRT precoder in (fTSl l. 
Then, we divide both the desired and the interference parts by 
A and 1/M^. The desired signal power is written as 


c-MRT 

^k,n 


1 H ~ 

j^2 Sk,n®k,n-^nSk,0 

1 

-^42{27rfuT,n)e-« +"?)" 


gk,0®k,nAlik,0 


( 22 ) 

(23) 


where is the mth diagonal element of the matrix 

expressing the variance of the mth element. Moreover, in the 
last step of (| 2 ^ . we have used the law of large numbers 
for large M in the case of collocated BS antennas, as well 
as that An is a scaled identity matrix or a scalar in the 
case of identical SLOs or CLO, respectively. In other words, 
gk,o has i.i.d. elements with variance As far as the 

interference is concerned, the first and third terms of dni) 
vanish to zero as M 00 by means of the same law. Thus, 
we have 


rMRT 

■^k,n 


M^PaX 


7^jg(27r/DTsn)e 






fc Jmm 


Mpd 


(24) 


where we have applied the law of large numbers to obtain 
A = from (O. Fi¬ 
nally, we substitute and the result in (1211 1 is 

obtained. ■ 

Remark 2: Clearly, (1211) shows that the achievable rate per 
user depends on the selection of q. Specifically, if q < 1/2, 
7 fc „ is unbounded. On the contrary, the SINR of user k reduces 
to zero, if q > 1/2. This clearly indicates that the transmit 
power of each user has been reduced over the required value. 
Most importantly, when q = 1/2, the SINR approaches a finite 
non-zero limit given by 


T E^En ^ 2/0 f rr, \ 
lk,n = - 4 —Jo(27r/D7sn)e 


-3 erf. 




[Rfc]mm- 


(25) 


As a result, both phase noise and user mobility reduce the 
effective SINR, but neither of them changes the power scaling 
law. 


IV. Numerical Results 

This section presents some representative numerical ex¬ 
amples illustrating the impact of channel aging due to the 
combined effects of phase noise and user mobility on the 
performance of very large MIMO systems with MRT. The 
correctness of the proposed results is validated by Monte-Carlo 
simulations. 

The simulation setup, based on long-term evolution (LTE) 
system specifications, considers a single cell with radius of 
R = 1000 meters, and it assumes a guard range of tq = 100 
meters specifying the distance between the nearest user and 
the BS. The BS, including M antennas, broadcasts to K 
users being uniformly distributed within the cell. While the 
theoretical analysis considers general R^, this section takes 
into account only for large-scale fading (lognormal shadowing 
and path-loss), i.e., R^ = 4 Im with 4 = Zkj {rkIroY, 
where Zk is a log-normal random variable with standard 
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Static environment 



Fig. 1. Simulated and deterministic equivalent downlink sum-rates with 
MRT precoder in a static environment with respect to different number of BS 
antennas, and various values of phase noise. 


deviation a (a = 8 dB) expressing the shadow fading effect, 
Tk denotes the distance between user k and the BS, and v 
(v = 3.8) is the path loss exponent. The training duration 
consists of T — K symbols, and the phase noises at the 
BS and user LOs are simulated as discrete Wiener processes 
by Q, with increment standard deviations in the interval 
0° — 2° II20I . The coherence time is Tc = 1/4/d = 1 ms, 
where /d = 250 Hz is the Doppler spread corresponding to a 
relative velocity of 135 km/h between the BS and the users, if 
the center frequency is assumed to be /c = 2 GHz. Moreover, 
given that the bandwidth for LTE-A is W = 20MHz, the 
symbol time is Tg = 1/(2W) = 0.025 fxs. 

One useful metric is the sum spectral efficiency denoted by 

K 

(26) 

fc=i 

In the following hgures, the “solid”, “dash”, and “dot” lines 
designate the analytical results with no phase noise, (j)rn,n = 
0°,‘Pk,n = 2°, and = Vk,n = 2°, respectively. The 

bullets represent the simulation results. 

Fig [U depicts the achievable sum-rate by varying the number 
of BS antennas for different values of phase noise in a static 
environment (u = 0 Km/h). By increasing the hardware 
imperfection in terms of phase noise, the sum-rate decreases. 
When the phase noise at the user side has a variance of 
= 2° the rate loss is 56% with the BS having M = 30 
antennas, but when M reaches 300, the loss is smaller, i.e., 
47%. An extra phase noise of 0° in the case of of CLO/ILOs 
setup leads to further losses of 50% and 47%, if M = 30 and 
M = 300, respectively. 

The variation of S with the normalized Doppler shift /dTs 
is presented in Fig. |2] when M = 60 for different values of 
phase noise starting from perfect LOs (no phase noise) to high 
phase noise (cr^^ = = 2°). The effect of Doppler shift 

dominates over the impact of phase noise. In fact, for low 
velocities in the order of 30 km/h equivalent to f^Ts ~ 0.2, 
the degradation due to phase noise starts to become insignih- 
cant, but then the achievable rate can become so low that it is 
inadequate for practical applications. In the same hgure, the 


M=60 



Fig. 2. Simulated and deterministic equivalent downlink sum-rates with 
MRT precoder for M = 60 with respect to the normalized Doppler shift, and 
various values of phase noise. 


Static environment 



Number of antennas M 

Fig. 3. Required transmit power to achieve 1 bit/s/Hz per user with MRT 
precoder in a static environment with respect to different number of BS 
antennas, and various values of phase noise. 


straight line, which is parallel to the horizontal axis, illustrates 
the sum-rate with imperfect CSI, but no channel aging (no 
user mobility or phase noise). Furthermore, as phase noise 
increases, the performance worsens. 

Fig. |3] illustrates how pd varies with the number of BS 
antennas in a static environment (/dT/ = 0 ), in order to 
achieve 1 bit/s/Hz. Specihcally, pd decreases considerably 
when we increase M. Especially, a closer observation shows 
a reduction in the transmit power by approximately 1.5dB 
after doubling the number of BS antennas, which agrees with 
previously known results e.g., m- The more severe is the 
phase noise, the higher the downlink power is. 
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V. Conclusions 

Channel aging, being a limiting factor in massive MIMO 
systems, considered initially only user mobility. However, 
phase noise, which innovates with time, is more severe in 
large antenna arrays than in convetional systems, since the 
former include low quality antenna elements. As a result, 
its contribution to channel aging is substantial. Hence, we 
obtained an interesting joint channel-phase noise model, and 
we derived the deterministic equivalent of the downlink SINR. 
Notably, we showed that user mobility affects more system’s 
performance than phase noise. Moreover, the phase noise 
plays a significant role on the achievable rate in low mobility 
conditions. Interestingly, massive MIMO are still favorable 
under generalized channel aging. 


Appendix A 
Proof of Theorem[2] 

The derivation of the deterministic SINR starts by deriving 
the deterministic equivalent of the normalization parameter A. 
Specifically, the normalization parameter can be written by 
means of ( fT3] ) as 

K 


A = 


E 


tr A„GoGgA„ 
K 


E 


Ei=l igA2gfc^ogfc,c 

K ^ M ” < 

. i=l 


(27) 


where we have applied ||4l Lem. 1]. Similarly, regarding the 
rest part of the desired signal power, we obtain 


—E 
M 




=—E 
M 

=—E 
M 

X—E 
M 


S/c ,0 0 


fl" ^n^k,n-^n^k 


w 


■trAlD 


k- 


(28) 


(29) 


The deterministic equivalent signal power Sk,n = 
limM->-oo Sk,n becomes by using (|2^ as 






(30) 


where 4 = tr A^D^. 

Now, we proceed with the derivation of each term of (fTTl i. 
Thus, the first term can be written as 


1 


var 


Artgfc.O 

1 


M 2 


E 


Angle, 0 


a.s. 


M—¥oo 


^0, (31) 


where the property of the variance operator var (x) = E[a;^] — 
E^[a;] together with ® have been used. Lemma 1 in H 
enables us to derive the deterministic equivalent as 


M 


var 


gfc,n®fc." Angfc.O 


where S[ = Atr A^Dfc (R^ - A^Dfc). 

Similarly, use of El Lem. 1] in the last term completes the 
proof. Specifically, we have 


M 2 


E 


gfc,n®fc,n Angj.C 


A 

M2 


tr0fe,nA„Dj0fc„Rfc, 


= Aa" 

M * ’ 


(33) 


since gk n and gi o are mutually independent. Note here that 
= Atr A2D,Rfe. 
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